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Abstract 

We discuss the nature of nonequilibrium phase transitions in the Hamiltonian Mean Field model 

^" 
O ■ using detailed numerical simulation of the Vlasov equation and molecular dynamics. Starting from 

fixed magnetization waterbag initial distributions and varying the energy, the states obtained after 

C^ ' a violent relaxation undergoes a phase transition from magnetized to non-magnetized states when 

-(— > ' 
c/3 ' 

going from lower to higher energies. The phase transitions are either first order or composed by a 

cascade of phase reentrances. This result is at variance with most previous results in the literature 



'^ ' mainly based in Lynden-Bell theory of violent relaxation. The latter is a rough approximation 

O ' and consequently not suited for an accurate description of nonequilibrium phase transition in long 

range interacting systems. 
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The physics of long range interacting systems is an active topic of investigation due to 

n 

the unusual and intriguing phenomenology they present [1]. A pair potential interaction is 
considered long range if it scales at greater distances as r~" with a < d, where r stands 
for the inter-particle distance and d the spatial dimension. This slow decaying interparticle 
potential is responsible for the coupling of distant components of the system, a condition not 
encountered in short range systems. A remarkable feature of these systems is that energy is 
non-additive and this opens up many nonintuitive phenomena, e. g. in the microcanonical 
ensemble it is possible to have negative specific heat and temperature jumps characteriz- 
ing first order phase transition. In this context canonical and microcanonical statistical 
ensembles can therefore be nonequivalent. Gravitational systems is another example that 
is largely studied 12 -pi in the microcanonical ensemble, and other systems not less impor- 

n Q 

tant that encompass different areas of physics, as plasmas [6| , wave-particle interactions [7|] 
and many others domains of application. A comprehensive review of the subject may be 
found in [l|. These systems also present uncommon dynamical features. Starting from an 
initial nonequilibrium configuration, these systems rapidly evolve by a violent relaxation 
to Quasi-Stationary States (QSS), where they stay trapped for long lasting times scaled 
as an increasing function of the number of constituent particles, and usually much longer 
than the time of observation that experimentalists are bound. Their structure was long 
ago recognized as non-Boltzmannian states, and are now properly interpreted in terms of 
stable steady states of the Vlasov equation and statistical equilibrium states in the sense of 



Lynden-Bell theory of violent relaxation 
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Recently, a number of researchers studied nonequilibrium phase transitions in the Hamil 



tonian Mean Field (HMF) model [19| in the context of Lynden-Bell theory |10|-|16[|. They 
consider initial waterbag states with a given magnetization and looked for the final magneti- 
zation after the violent relaxation. They then observed a phase transition from a magnetized 
to a non- magnetized QSS. Nevertheless the nature of such phase transitions and whether 
Lynden-Bell theory correctly predicts them is still open to debate 12, 113]. 

In this paper we attempt to provide a more detailed description of the nature of nonequi- 
librium phase transitions of the HMF model, and in particular, we pay attention to reentrant 
phases that seem to play an important and previously not fully acknowledged role. We pro- 
vide results from numerical simulations of the Vlasov Equation and Molecular Dynamics 
(MD). The HMF model is a system of identical particles on a circle with unit mass and 



Hainiltonian: 
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where ^j is the angle that particle i makes with a reference axis and Pi stands for its conjugate 
momentum. The 1/A^ factor in the potential energy corresponds to the Kac prescription to 
make the energy extensive and justify the validity of the mean field approximation in the 
limit A^ — ?► oo. The relevant order parameter is the magnetization defined as: 



M = JM^ + M^, 



(2) 



where M^ = (l/N) E» cos 9i and My = (1/iV) J2i sin Oi. 

In the continuum limit the evolution of the single particle distribution function f{6,p,t) 
is governed by the Vlasov equation |l[ |l8| : 

df df _dV[f]df _ 

where V[f] is is the interaction potential that depends self-consistently on f{9,p,t) and is 
given by V[f] (6) = I - M^f] cos (6) - My[f] sin (6), with 



M,[f]= dedpf{9,p,t)cos{9) 



and 



My[f]= d9dpf{9,p,t)sm{9). 



(4) 



(5) 



In the foregoing discussion and following previous approaches 10Nl6| we consider as initial 
state a waterbag distribution, i. e. f{9,p,t = 0) = l/2ApA9 if < 6* < A9 and \p\ < Ap, 
and f{9,p,t = 0) = otherwise. The initial magnetization and energy (per particle) are 
given by M = [(1 - cos A^)^ + (sin A9f] ^^^ /A9, and e = ApV24 + (1 - M^)/2. 

In order to discuss the out of equilibrium phases corresponding to the final state after 
a violent relaxation, it is important to establish how long it takes for the system to settle 
down into a QSS or a possibly perpetually oscillating steady state [20|. Figure [1] shows 
the magnetization as a function of time for different initial values of magnetization Mq and 
energy e. It becomes clear that in many cases the QSS in only attained (or approached) for 
times of order 10'^, at least one order of ma gnit ude greater than the total time used in some 
previous simulations on the same problem 2ll-l23|. 



Antoniazzi et al. 



21| compared the predictions from Lynden-Bell theory with N-body 



simulations at ii^ = 0.69 and obtained a reasonable agreement for Mq < 0.897. They also 




FIG. 1. Mono-Log graph of magnetization as a function of time for some initial magnetizations 
and energies per particle. 



constructed a phase diagram in the (Mq, e) plane and noticed that the system presents 
first and second order phase transitions separated by a tricritical point. They proceed 
further and performed numerical simulations of the Vlasov equation and found reasonable 
agreement with N-body simulation and Lynden-Bell theory. Staniscia et al. 22[ in their 
calculations confirmed the existence of reentrant phases as predicted by theory but show 
anyhow some discrepancies and argued that this occurs due to incomplete relaxation during 
violent relaxation. 

Figures [2H1] show the final magnetization as a function of energies for a few representa- 
tive values of the initial magnetization Mq from the solution of the Vlasov equation, MD 



simulations and Lynden-Bell theory 2J|. Vlasov simulations were performed using a Vlasov 



integrator code in Ref. 



with a numeric grid with 512 x 512 points in the one particle 
phase space, total integration time t/ = 3000.0 and averaging from t = 2000.0 to t = 3000.0. 
For the more detailed graphics in figures [2b, EJi, [3b, [lb and [5b we used a 2048 x 2048 grid 
with integration time tf = 1000.0 and averaging from t = 800.0 to tf. The results from 
Lynden-Bell theory were obtained using the approach in 26| . Figures [2^ and ^ indicate 
that the transition is discontinuous in both cases predicted be from Lynde-Bell theory for 
Mo = 0.1 but not for Mq = 0.3. This is even more clearly shown in Figs. [2b and [2jl that 
show the region near the phase transition using more simulation points and more accurate 
Vlasov simulations. For Mq = 0.1 at least three reentrant phase transitions are observed 
before the predicted (and observed) phase transition. The phase transition is more neatly 
observed for Mq = 0.3, where no mater the order parameter chosen (M or M^) the transition 



is clearly first order, from both MD and Vlasov equation solution. In fact some discussion 
exists in the literature whether M or M^ should be used as an order parameter [lO(] . Here we 
argue that both choices lead to the same characterization of the order of the phase transi- 
tions. For Mq = 0.4 the discontinuity in the phase transition is even more evident as shown 
if Fig. El We note that for this particular value of magnetization the discontinuity in the 



phase transition was previously reported by Pakter and Levin [17|. They were also able to 
correctly predict the phase transition using a new ansatz for the distribution function based 
on dynamical properties of the underlying Hamiltonian dynamics. The situation gets even 
more interesting for Mq = 0.5 in Fig.HJwhere a close look around the phase transition reveals 
a cascade of reentrant phases. As a consequence it is not clear how to asses the nature of the 
phase transition in this case. For higher initial magnetizations the same analogous behavior 
is observed. 

It is important to note that the critical energy of the phase transitions as predicted by 
Lynden-Bell theory is only an approximation, albeit a good one. On the other hand, all 
previous studies of nonequilibrium phase transitions in the HMF model have concentrate 
on the magnetization as an order parameter, which is obtained from the spatial distribution 
function. It is interesting also to discuss what occurs with the velocity distribution function 
along the same lines depicted previously. For that purpose we use the moments of the 
velocity distribution function given by the average of powers of v a.s fik = {v'^)- Figure [5] 
shows the averaged moments /i4 and /le for Mq = 0.4. It is quite evident that Lynden-Bell 
theory gives reasonable results only for lower energies. The right panel of the same figures 
shows a discontinuity in fi^ and /le, a clear indication that the phase transition is indeed 
first order. 

In this paper we have investigated phase transitions of QSS's using three different ap- 
proaches: Lynden-Bell theory of violent relaxation, numeric solutions of the Vlasov equation 
and molecular dynamics. Previous points in favor of Lynden-Bell theory is that it gives a 
reasonable first approximation of the QSS's, and in this context, also allows to predict out- 
of-equilibrium phase transitions, although it is also accepted the argument that the QSS's 
are, or can be, incomplete mixed stable states of the Vlasov equation J27|. The results 
presented here show unequivocally that the nature of phase transitions is of first order are 
noticeable for different magnetizations, and reentrant phases are more common than pre- 
viously noted, as the cascade of phase reentrances observed for Mq = 0.1 and Mq = 0.5 
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FIG. 2. QSS magnetization Mj as a function of the energy per particle for initial magnetizations 
Mq = 0.1 (panels a,b) and Mq = 0.3 (panels (c,d). All points from Vlasov solution were obtained 
using a numeric grid with 512 x 512 points and a time step At = 0.2, total integration time 
t/ = 3000.0, and averaging from t = 2000.0 up to t = tf except (b) that used a grid with 
2048 X 2048 points, tf = 1000.0 and averaging from t = 800.0 up to tf. Molecular Dynamics 
simulations (MD) were performed with N = 20, 000, 000. 

clearly illustrates. The simulations also show that Lynden-Bell theory is not suitable to 
accurately predict these transitions. Molecular dynamics results are in very good agreement 
with numeric solutions of the Vlasov equation. As a step forward we have decided to calcu- 
late the moments (4th and 6th) of the velocity distribution and once more they diverge of 
those predicted by Lynden-Bell theory. Although the latter, according to our calculations, 
is inadequate to explain nonequilibrium phase transitions in QSS's it predicts with some 
accuracy the position of the phase transition, but not its order, and certainly not the phase 
reentrances here reported. At lower energies it yields quite reasonable results for magne- 
tization and moments of the velocity distribution function, but strongly depart from the 
correct values at higher energies. Therefore a detailed and accurate study of nonequilibrium 
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FIG. 3. QSS magnetization for Mq = 0.4 computed from Lynden-Bell theory (LB), numeric 
solution of Vlasov equation and Molecular Dynamics (MD) with N = 20, 000, 000, t/ = 3000.0. 
The right panel shows in greater detail the discontinuity in the magnetization. 
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FIG. 4. Final magnetization as a function of energy for Mq = 0.5 from the solution of Vlasov 
equation. The left panel (a) also shows the prediction from Lynde-Bell theory (LB). 



phase transitions in long-range interacting systems cannot be based on Lynden-Bell theory. 
Unfortunately a completely satisfactory theory for violent relaxation is still lacking, even 
though some progress was obtained in Refs. 28| and 29 |. 
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FIG. 5. Fourth and sixth moments of the velocity distribution function for Mq = 0.4 as a function of 
energy computed from Lynde-Bell theory (LB), the solution of the Vlasov equation and Molecular 
Dynamics (MD) with N = 20, 000, 000. 
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